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Abstract 

The lifetimes of B + and B° mesons have been measured using the entire 
sample of 550,000 hadronic Z° decays collected by the SLD experiment at 
the SLC between 1993 and 1998. In this paper, we describe the inclusive 
analysis of the 350,000 hadronic Z° decays collected in 1997-98 with the up- 
graded SLD vertex detector. In this data period, a high statistics sample 
of 30903 (20731) charged (neutral) vertices with good charge purity is ob- 
tained. The charge purity is enhanced by using the vertex mass, the SLC 
electron beam polarization (73% for 1997-8) and an opposite hemisphere jet 
charge technique. Combining the results of this data sample with the re- 
sults from the earlier data yield the following preliminary values: t b + = 
1.623 ± 0.020(stat) ± 0.034(syst) ps, t b o = 1.589 ± 0.021(stat) ± 0.043(syst) 
ps, t b +/t b o = 1.03T ±8:811 (stat) ±0.024(syst). 
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The spectator model predicts that the lifetime of a heavy hadron depends 
upon the properties of the constituent weakly decaying heavy quark Q and 
is independent of the remaining, or spectator, quarks in the hadron. This 
model fails for the charm hadron system for which the lifetime hierarchy 
r D + ~ 2.?>r D + ~ 2.5t d o ~ 5r A + is observed. Since corrections to the spectator 
model are predicted to scale with l/rug the B meson lifetimes are expected to 
differ by less than 10% Q. Hence a measurement of the B + and -B lifetimes 
provides a test of this prediction. In addition, specific B meson lifetimes are 
needed for many important measurements, e.g. to determine the element V c b 
of the CKM matrix. 

An analysis has been reported using the 1993-5 data sample of 150,000 
hadronic Z° decays collected with the original CCD vertex detector (VXD2) 
as well as with the first 50,000 hadronic Z° decays collected in 1996 using the 
upgraded vertex detector ( VXD3) M . In this paper we describe the analysis 
of the additional 350,000 hadronic Z° decays collected in the 1997-1998 run 
by the SLD detector at the SLC, and combine these results with those from 
the previous work||. The excellent 3-D vertexing capabilities of SLD are 
exploited with an inclusive topological vertexing technique Q to identify B 
hadron vertices produced in hadronic Z° decays with high efficiency (This 
inclusive technique has the advantage of very efficient B vertex reconstruction 
since most B decays are used). The decay length is measured using the 
reconstructed vertex location while the B hadron charge is determined from 
the total charge of the tracks associated with the vertex. Knowledge of the 
average energy of 5-hadrons produced in Z° decays, in conjuction with the 
decay length, lets one infer the B lifetime. 

The components of the SLD utilized by this analysis are the Central 
Drift Chamber (CDC)|| for charged track reconstruction and momentum 
measurement and the CCD pixel Vertex Detector (VXD)0, ^[ for precise po- 
sition measurements near the interaction point. These systems are immersed 
in the 0.6 T field of the SLD solenoid. Charged tracks reconstructed in the 
CDC are linked with pixel clusters in the VXD by extrapolating each track 
and selecting the best set of associated clusters 0. For VXD3 the track im- 
pact parameter resolutions at high momentum are 9 /xm and 11 jum in the 
rcf) and rz projections respectively (z points along the beam direction), while 
multiple scattering contributions are 33 /xm /(psin 3//2 #) in both projections 
(where the momentum p is expressed in GeV/c). 

The centroid of the micron-sized SLC Interaction Point (IP) in the rip 
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plane is reconstructed with a measured precision of aip = (4 ± 2) /iin using 
tracks in sets of ~ 30 sequential hadronic Z° decays. The median z position 
of tracks at their point of closest approach to the IP in the rip plane is used 
to determine the z position of the Z° primary vertex on an event-by-event 
basis. A precision of ~ 30 /mi on this quantity is estimated using Z° — > bb 
Monte Carlo simulation for VXD3. 

The simulated Z° — > qq events are generated using JETSET 7.4 ||. The 
B meson decays are simulated using the CLEO B decay model tuned to 
reproduce the spectra and multiplicities of charmed hadrons, pions, kaons, 
protons and leptons as measured at the T(4S) by ARGUS and CLEO [0]. 
The branching fractions of the charm hadrons are tuned to the existing mea- 
surements [[5J. The B mesons and baryons are generated with lifetimes of 
t b + = 1.64 ps, t b o = 1.55 ps, t b o = 1.57 ps, and T\ b = 1.22 ps. The 6-quark 
fragmentation follows the Peterson et al. parameterization |TIJ; the mean 
value of the fragmentation function in the MC generation was 0.698. The 
SLD detector is simulated using GEANT 3.21 

Hadronic Z° event selection requires at least 7 CDC tracks which pass 
within 5 cm of the IP in z at the point of closest approach to the beam and 
which have momentum transverse to the beam direction p? >200 MeV/c. 
The sum of the energy of the charged tracks passing these cuts must be 
greater than 18 GeV. These requirements remove background from Z° — > 
l + l~ events and two-photon interactions. In addition, the thrust axis de- 
termined from energy clusters in the calorimeter must have |cos 6*| < 0.85. 
within the acceptance of the vertex detector. These requirements yield a 
sample of ~ 267, 500 hadronic Z° decays for the 1997-98 datset. 

Good quality tracks used for vertex finding must have a CDC hit at a 
radius<39 cm, and have >23 hits to insure that the lever arm provided by 
the CDC is appreciable. The CDC tracks must have px >250 MeV/c and 
extrapolate to within 1 cm of the IP in r<p and within 1.5 cm in z to eliminate 
tracks which arise from interaction with the detector material. The fit of the 
track must satisfy x 2 /d.o.f.< 8. At least two good VXD3 links are required, 
and the combined CDC/VXD fit must also satisfy x 2 /d.o.f.< 8. 

The topological vertex reconstruction is applied separately to the tracks 
in each hemisphere (defined with respect to the event thrust axis). The 
vertexing algorithm is described in detail in Ref. and summarized here. 
The vertices are reconstructed in 3-D coordinate space by defining a vertex 
function V(r) at each position r. The helix parameters for each track i are 
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used to describe the 3-D track trajectory as a Gaussian tube fi(r), where the 
width of the tube is the uncertainty in the measured track location close to 
the IP. A function f (r) is used to describe the location and uncertainty of 
the IP. V(r) is defined as a function of /o( r ) an d the /j(r) such that it is small 
in regions where fewer than two tracks (required for a vertex) have significant 
/i(r), and large in regions of high track multiplicity. Maxima are found in 
V(r) and clustered into resolved spatial regions. Tracks are associated with 
these regions to form a set of topological vertices. 

The efficiency for reconstructing at least one secondary vertex in a b hemi- 
sphere is ~ 67% for VXD3. The efficiency falls at shorter decay length as it 
becomes harder to resolve the secondary vertex from the IP. For hemispheres 
containing secondary vertices, the 'seed' vertex is chosen to be the one with 
the highest V(r) value. Vertices consistent with a If ° — > it + -k~ decay, within 
a 14 MeV invariant mass window, are excluded from the seed vertex selection 
and the two tracks are discarded. 

A vertex axis is formed by a straight line joining the IP to the seed 
vertex. The 3-D distance of closest approach of a track to the vertex axis, 
T, and the distance from the IP along the vertex axis to this point, L, are 
calculated for all quality tracks. Monte Carlo studies show that tracks which 
are not directly associated with the seed vertex but which pass T< 0.1 cm 
and L/D> 0.3 (where D is the distance from the IP to the seed vertex) are 
more likely to have been produced by the B decay sequence than to have an 
alternative origin. Hence such tracks are added to the set of tracks in the seed 
vertex to form the candidate B decay vertex, containing tracks from both 
the B and cascade D decays. This set of tracks is fitted to a common vertex. 
If the probability of the vertex fit is greater than 5% the B decay location is 
taken to be the fitted vertex location. Since the vertex includes tracks form 
both the B and cascade charm decay points the fit probability distribution 
is not flat and ~ 65% of vertices reconstructed have a fit probability below 
5%. The reconstruction of the B decay location in this case is improved 
by dividing the tracks (if there are at least three) into two 'sub-vertices' 
by selecting the combination of tracks with maximum two sub-vertex fit 
probability. The B decay location is given by the sub-vertex closest to the 
IP, reducing systematic dependence upon the physics of the B — > D decay. 
The distance from the IP to the B decay location is the reconstructed decay 
length. Since the purity of the B charge reconstruction is lower for decays 
close to the IP, where tracks are more likely to be wrongly assigned, decay 
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lengths are required to be > 1 mm. To avoid using the vertices of tracks 
originating from interactions with the detector material, the distance of the 
decay vertex from the beamline is required to be < 22 mm, i.e. more than 
1 mm inside the SLC beampipe. 
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Figure 1: (a) Mp T of reconstructed vertex for 1997-98 data (points) and 
Monte Carlo (histogram), (b) B + fraction in the charged sample, (c) B° 
fraction in the neutral sample. 

The mass M of the reconstructed vertex is calculated by assuming each 
track has the mass of a pion. The transverse component Pt of the total 
momentum of vertex tracks relative to the vertex axis is calculated in order 
to determine the Pt corrected mass: 
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Mp T = y M 2 + Pj, + | Pj- 1 . (1) 

This quantity is the minimum mass the decaying hadron could have in 
order to produce a vertex with the quantities M and Pp. The direction of 
the vertex axis is varied within the la limits constraining the axis at the 
measured IP and reconstructed seed vertex such that the Pt is minimized 
within this variation. This procedure prevents non-S background vertices 
acquiring a high Mp T due to a fluctuation in the measured Pt- The accurate 
3-D vertexing and precisely measured IP at SLD allow significant gain in the 
6-tag efficiency with high purity using this technique JE2 . 



A comparison of the distribution of Mp T in data and Monte Carlo is 
shown in Fig. |](a). Vertices from decays surviving the rejection can 
be seen around 0.5 GeV/c 2 . This figure shows that a large fraction of the 
charm and light flavor contamination in the sample is eliminated by requiring 
Mp T > 2 GeV/c 2 . It is also required that Mp T < 5.2 GeV/c 2 since vertices 
with Mp T greater than the B meson mass are likely to contain background 
tracks. These cuts yield a sample with b hemisphere purity of ~ 98% with 
an efficiency of 38%. The average B decay vertex multiplicity is 5.0 tracks. 

To improve the B hadron charge reconstruction, tracks which fail the 



initial selection but have pp > 200 MeV/c and yer^ + o 2 z < 700 /im, where 
c r <p {crz) is the uncertainty in the track position in the rip (rz) plane close to 
the IP, are considered as decay track candidates. The charge of these tracks 
which pass the cuts T< 0.1 cm and L/D> 0.3 is added to the B decay 
charge. On average, 0.4 tracks pass these criteria in b hemispheres. These 
lower quality tracks are used only to improve the charge reconstruction. 

Fig. ^ shows a comparison of the reconstructed charge between data and 
Monte Carlo for the 1997-1998 dataset. At this stage the charged sample 
consists of 30903 vertices with vertex charge equal to ± 1,2 or 3, while the 
neutral sample consists of 20731 vertices with charge equal to 0. Monte 
Carlo studies indicate that the charged sample is 97.0% pure in B hadrons 
consisting of 57.2% B+, 32.0% B°, 8.1% B°, and 4.5% B baryons. (Charge 
conjugation is implied throughout this paper with the exception of the no- 
tation B + and _B_ introduced later to distinguish the charged mesons B(bu) 
and B(bu) and respectively.) Similarly, the neutral sample is 98.3% pure in 
B hadrons consisting of 21.9% B+, 55.0% B°, 15.6% B° s and 7.6% B baryons. 
The statistical precision of the measurement depends on the separation be- 
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Figure 2: Reconstructed vertex charge for data (points) and Monte Carlo 
(histogram) for the 1997-98 dataset. 

tween the B + and B° in these samples. 

The lifetime measurement relies on the ability to separate B + and B° 
decays by making use of the vertex charge. Monte Carlo studies show that the 
purity of the charge reconstruction is more likely to be eroded by losing tracks 
from the B decay chain through track selection inefficiencies and track mis- 
assignment than by gaining mis-assigned tracks originating from the primary 
or other background to the B decay. Furthermore, the decays which are 
missing some B tracks tend to have lower vertex mass as well as lower charge 
purity. Fig. |l|(b) and (c) show the fraction of B + decays in the charged sample 
and the fraction of B° decays in the neutral sample respectively for MC. The 
probability for the vertex to originate from a positively charged, neutral 
or negatively charged B meson is denoted as P C (B + ), P C (B°) and P C (BJ) 
respectively, normalized such that P C (B + ) + P C (B°) + P C (B_) = 1. Vertices 
are given a weight, w, according to their analyzing power for separating B + 
and B° decays: w = \2P C (B + ) - 1|, where P C (B+) = P C (B+) + P c (5__). The 
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probabilities, and hence the weight, are a function of Mp T . 

The charge reconstruction in Fig. ^| shows good agreement between data 
and MC. A further check is made using the SLC electron beam polarization. 
The polarized forward-backward asymmetry Aps{Pe, cos#) can be described 
by 

^ ( P e , C039) . 2A ^_^_, (2) 

where Ab = 0.94 and A e = 0.155 (Standard Model values), P e is the elec- 
tron beam longitudinal polarization, and 9 is the angle between the thrust 
axis and the electron beam direction (the thrust axis is signed such that it 
points in the same hemisphere as the reconstructed vertex). Using negative 
(positive) vertex charge, with vertices weighted by the Mp T dependent an- 
alyzing power, to tag the b (b) quark flavor the resulting forward backward 
asymmetry is sensitive to the accuracy of the vertex charge reconstruction. 
Good agreement between data and MC can be seen in Fig. |] for the 1997-8 
data, indicating that the MC adequately reproduces the charge reconstruc- 
tion purity of the data. (Random vertex charge assignment would generate 
distributions with no signed cos# asymmetry.) 

The polarized forward-backward asymmetry can be used to tag the initial 
state b or b flavor of the hemisphere. The initial state b/b tag is also used to 
enhance the charged sample purity by giving a higher (lower) weight to the 
B + hypothesis if the vertex charge agrees (disagrees) with the b/b tag. The 
probability for correctly tagging a b quark at production using the e~ beam 
polarization is expressed as 

Pa® = 1 + A ^ cos ^ . (3) 

A jet charge technique is used in addition to the polarized forward-backward 
asymmetry. For this tag, tracks in the hemisphere opposite that of the recon- 
structed vertex are selected. These tracks are required to have momentum 
transverse to the beam axis p± > 150 MeV/c, total momentum p < 50 
GeV/c, impact parameter in the plane perpendicular to the beam axis 5 < 2 
cm, distance between the primary vertex and the track at the point of closest 
approach along the beam axis Az < 10 cm, and | cos#| < 0.87. With these 
tracks, an opposite hemisphere momentum-weighted track charge is defined 
as 
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Figure 3: Distribution of cos# signed by (P e x vertex charge) for 1997-8 data 
(points) and Monte Carlo (histogram). 
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where % is the electric charge of track i, its momentum vector, f is the 
thrust axis direction, and k is a coefficient chosen to be 0.5 to maximize the 
separation between b and b quarks. The probability for correctly tagging a b 
quark in the initial state of the vertex hemisphere can be parameterized as 



P Q (b) 



(5) 



2_ — |— g&Qopp ' 

where the coefficient a = —0.27 as determined using the Monte Carlo sim- 
ulation. This technique is independent of the polarized forward-backward 
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asymmetry tag. The average purity of the b/b tag is ~ 72% using the forward 
backward asymmetry (\P e \ = 73%) and 67% for the jet charge technique. 

The two initial state tags can be combined to form an overall initial state 
tag with b/b quark probability Pj(6)/Pj(6) (a function of P e , cos 6* and Q opp ). 
This probability is then combined with that obtained from the vertex charge 
reconstruction (a function of Mp T ) to determine the overall probability of a 
B + or B° decay: 

p(B+) = P c (B + )P(b) + P c (B_)P t (b) 

1 ' P c (B + )P(b) + P c (B.)P(b) + 0.5 x P C (B°) { ' 

where B + and B_ denote the positively and negatively charged B mesons 
separately and P(B°) = 1 - P(B+). If the probability P(B+) > P(B°) the 
vertex is classified as charged, otherwise it is added to the neutral sample. 
In either case it is weighted by the analyzing power w = \2P(B + ) — 1|. 
Including the initial state tag information enhances the statistical power of 
the analysis by ~20%. 
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Figure 4: Decay length distributions for the 1997-8 data (points) and best 
fit Monte Carlo (histogram). 

The B + and B° lifetimes are extracted from the decay length distributions 
of the B vertices in the charged and neutral samples using a binned x 2 
fit. These distributions are fitted simultaneously to determine the B + and 
B° lifetimes. For each set of parameter values, Monte Carlo decay length 



10 



distributions are obtained by reweighting entries from generated B + and B° 
decays in the original Monte Carlo decay length distributions with W(t, r) = 
(~ e~*/ T ^) / (t^ - e - */ 7 " 96 "), where r is the desired B + or B° lifetime, r gen is 
the average lifetime generated in the Monte Carlo, and t is the proper time 
of each decay The fit then compares the decay length distributions from 
the data with the reweighted Monte Carlo distributions. Fig. |] shows the 
reconstructed decay length for data and best fit MC for the 1997-8 data and 
the charged and neutral samples. (The fits are made using histograms in 
which the bin size increases with decay length such that the number of entries 
per bin remains approximately constant.) The fit to the 1997-8 data sample 
yields lifetimes of r B + = 1.657 ± 0.023 ps and t b o = 1.578 ± 0.024 ps, with a 
ratio of t b +/t b o = 1.049±g;g|f and a xV d -°- f - = 79.8/76. These results are 
then corrected to correspond to b-fragmentation with a mean scaled energy 
< X B >= 0.714 ||13|| , a B® lifetime of 1.49 ps, and a b-baryon fraction of 
10.2% @. The corrected results for 1997-98 are: t b + = 1.613 ±0.023 ps and 
t b o = 1.565 ± 0.024 ps, with a ratio of t b +/t b o = 1.030^7- 

To account for a discrepancy between data and Monte Carlo in the frac- 
tion of tracks passing the selection criteria, a ~4% tracking efficiency correc- 
tion with dependence on track momenta and angles is applied to the simu- 
lation [EN. The corrected Monte Carlo is used in the lifetime fits, with the 
effect of the entire correction taken as the systematic error. 

The physics modeling systematic uncertainties were determined as fol- 



lows. The mean fragmentation energy <xe> of the B hadron [13] and 



the shape of the x B distribution [|TJJ] were varied. Since the fragmenta- 
tion is assumed to be identical for the B + and B° mesons, this uncer- 
tainty has little effect on the lifetime ratio. The four branching fractions for 
B + /B° — > D°/D~X were varied by twice the uncertainty given in Ref. [15 



for B — > D°/D X. The fraction of B + /B° decays producing a DD pair 
was also varied. The average B + and B° decay multiplicity was varied by 
±0.3 tracks |6 in an anticorrelated manner. Uncertainties in the B® and B 
baryon lifetimes and production fractions mostly affect the B° lifetime since 
the neutral B° s and B baryon are a more significant background for the B° 
decays. The systematic errors due to uncertainties in charmed meson decay 
topology were estimated by changing the Monte Carlo D decay charged mul- 
tiplicity and K° production according to the uncertainties in experimental 
measurements ||17|| . The effect of varying the lifetime of charm hadrons {D + , 
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D°, D s , A c ), as well as their momentum spectra in the B decay rest frame 
was found to be negligible. 

The fitting uncertainties were determined by varying the bin size used in 
the decay length distributions, and by modifying the cuts on the minimum 
decay length (no cut-2 mm) and maximum radius cuts (20 mm-no cut) used 
in the fit. Fit results are consistent within statistics for these variations, but 
a systematic error is conservatively assigned using the RMS variation of the 
results. 

Table [l] summarizes the systematic errors on the B + and B° lifetimes and 
their ratio for the 1997-1998 dataset. 

Finally, the 1993-95 and 1996 results || are then combined with these 
more recent results, taking into account correlated uncertainties. For the 
earlier data sets, the results have been adjusted so that their assumed values 
for the b-baryon fraction and the mean of the fragmentation function agrees 
with this current analysis. The combination of the 1993-5, 1996, and 1997-98 
data samples yield lifetimes of r B + = 1.623±0.020 ps and t b o = 1.585 ±0. 031 
ps, with a ratio of t b +/t b o = 1.037t K o24- The combined systematics for the 
1993-98 data are shown in Table [| 

In summary, from the entire 550,000 Z° decays collected by SLD between 
1993 and 1998, the B + and B° lifetimes have been measured using an inclu- 
sive topological technique. The analysis of the 1997-98 dataset of 350,000 de- 
cays, discussed above, isolates 51634 B hadron candidates with good charge 
purity enhanced by the vertex mass, e~ beam polarization and opposite hemi- 
sphere jet charge information. These results have been combined with the 
earlier 1993-96 measurements to yield the following preliminary results: 



r B+ = 1.623 ± 0.020(stat) ± 0.034(syst) ps, 
r B o = 1.585 ± 0.021(stat) ± 0.043(syst) ps, 

— = 1.037 ±2 024 (stat) ±0.024(syst). 

T B 

These results are consistent with the expectation that the B + lifetime is 
up to 10% greater than that of the B° and have the best statistical accuracy 
among current measurements 
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Table 1: Summary of systematic uncertainties in the B + and B° lifetimes 
and their ratio for the 1997-98 dataset. 



Systematic Error 




Ar B+ 


At b o 


T B 






(ps) 


(ps) 


Detector Modeling 


Tracking efficiency 




0.005 


0.011 


0.008 


Tracking resolution 




0.003 


0.003 


0.003 




Physics 


Modeling 






b fragmentation 


0.714 ± 0.008 


0.025 


0.030 


<0.003 




xe shape 


0.011 


0.009 


<0.003 


BR(B -> DX) 




0.005 


0.009 


0.007 


BR(5 -> DDX) 


0.19 ± 0.05 


0.015 


0.014 


0.016 


B decay multiplicity 


5.3 ± 0.3 


<0.003 


<0.003 


<0.003 


B° s fraction 


0.115 ± 0.020 


0.003 


<0.003 


<0.003 


B baryon fraction 


0.102 ± 0.020 


0.004 


0.017 


0.008 


Bg lifetime 


1.49 ± 0.06 ps 


<0.003 


0.014 


0.009 


B baryon lifetime 


1.22 ± 0.05 ps 


<0.003 


0.005 


0.003 


D decay multiplicity 




0.003 


0.005 


0.006 


D decay K° yield 




<0.003 


0.008 


0.004 




Monte Carlo and Fitt 


ing 




Fitting systematics 




0.010 


0.009 


0.006 


MC statistics 




0.006 


0.006 


0.007 


TOTAL 




0.033 


0.044 


0.026 
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Table 2: Summary of systematic uncertainties in the B + and B° lifetimes 
and their ratio for the combined 1993-8 data. 



Systematic Error 




Ar B+ 


At b o 


T B 






(ps) 


(ps) 


Detector Modeling 


Tracking efficiency 




0.004 


0.004 


0.005 


Tracking resolution 




0.003 


0.003 


0.003 




Physics 


Modeling 






b fragmentation 


0.714 ± 0.008 


0.025 


0.028 


.004 




xe shape 


0.011 


0.010 


<0.003 


BR(B -> DX) 




0.005 


0.009 


0.007 


BR(5 -> DDX) 


0.19 ± 0.05 


0.013 


0.012 


0.014 


B decay multiplicity 


5.3 ± 0.3 


<0.003 


<0.003 


<0.003 


B° s fraction 


0.115 ± 0.020 


0.005 


<0.003 


<0.003 


B baryon fraction 


0.102 ± 0.020 


0.005 


0.017 


0.008 


Bg lifetime 


1.49 ± 0.06 ps 


<0.003 


0.014 


0.009 


B baryon lifetime 


1.22 ± 0.05 ps 


<0.003 


0.005 


0.003 


D decay multiplicity 




0.005 


0.006 


0.008 


D decay K° yield 




0.003 


0.009 


0.006 




Monte Carlo and Fitt 


ing 




Fitting systematics 




0.008 


0.007 


0.005 


MC statistics 




0.005 


0.005 


0.006 


TOTAL 




0.034 


0.043 


0.024 
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